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ABSTRACT 

- Method and device for neutralizing, by controlled gas injection, the formation of 
liquid slugs at the foot of a pipe portion greatly inclined to the horizontal or riser 
connected to a pipe canying circulating multiphase fluids such as hydrocarbons. 

- Flow^ rate control is essentially obtained by injecting, at the base of the riser, a 
volume of gas proportional to the flow rate variation with time of the gas phase of the 
circulating fluids, and preferably substantially equal thereto, when this variation is 
positive. This action can be completed by modulating the injected gas flow by a 
quantity proportional to the flow rate variation of the liquid phase of the circxxlating 
fluids. 

- Applications : offshore hydrocarbon production control for example. 
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FIELD OF THE INVENTION / 

The present invention relates to a method and to a device for neutralizmg, by 
controlled gas injection, liquid slugs or accumulations at the foot of a pipe portion 
greatly inclined to the horizontal or riser connected to a pipe carrying circulating 
multiphase fluids such as hydrocarbons. 

BACKGROUND OF THE INVENTION 

In order to make deep-sea reservoirs or marginal fields sufficiently cost-effective, 
oil companies have to develop new development techniques, as econoroical as possible. 
It is thus more advantageous to directly transport the two-phase mixture consisting of 
liquid (oil and a httle water) and gas in a single pipe or pipeline to onshore fecilities in 
order to be separated. A pipe portion greatly inclined to the horizontal (often close to 
the vertical), referred to as riser by specialists, which is connected to the deep-sea pipe, 
is used therefore. However, the gas and the liquid being transported together, flow 
instability phenomena may occur in the zone of connection with tihie riser, which lead to 
serious development problems. 

In particular, when the gas and liquid inflow rates are low, the liquid phase 
accumulates in the lower parts of the pipeline and stops the gas flowing past. The 
upstream pressure increases and eventually expels the liquid slug to another low part or 
even in the phase separator at the outlet. These accumulation phenomena can reduce the 
productivity and fill pipes designed to receive gas with liquid. One of these phenomena, 
commonly referred to as severe slugging by specialists, has formed the subject of many 
studies, either experimental by means of test loops, or by simulation by means of 
simulation softwares such as, for example, the TACITE simulation code which is 
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notably the object of the following patents or patent applications : US~5,550,76l5 FR- 
2,756,044 (US-6,028,992) and FR-2,756,045 (US-5,960,187X EEl-00/08,200 and FR- 
00/09,889 fUed by the applicant. 

This slugging phenomenon is described hereafter in the simple case illustrated in 
5 the accompanying figures where a pipe of low inclination and a riser ended by a 
separator designed to separate the liquid phase jfrom the gas phase are considered. 

The liquid accumulates in the lower part of the pipe and tends to stop the gas 
flowing past. The gas is conq>ressed until the upstream pressure exceeds the pressure 
due to the weight of the accumulated liquid. A long liquid slug is then pushed by the 
10 expanding gas. Under such conditions, an alternating phenomenon is observed, where 
the liquid blocks tiie gas phase, then flows off under the pressure of the gas and 
eventually accumulates and blocks the gas again. 

More precisely, the periodic process takes place as follows : 

Stage I : the liquid accumulates at the foot of the riser and stops the gas flowing 
15 past. The pressure rises ; 

Stage II : the upper level of the liquid having reached the top of the riser, the liquid 
phase flows into the separator ; 

Stage III : the gas pocket reaches the foot of the riser and flows into the riser. The 
slug flows into the separator with a much higher velocity ; the gas pocket « explodes » 
20 in the riser ; 
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Stage IV : when the gas pocket reaches the top of the riser, the pressure at the foot 
of the pipe has a minimum value. The liquid falls down along the wall of the riser. It 
accumulates again at the foot of the riser and a new cycle starts. 

A well-known technique referred to as gas lift by specialists allows to overcome 
5 this phenomenoa It essentially consists in permanently injecting gas at the base of the 
riser to prevent the accumulation of liquid at the bottom of the pipe. Since this 
phenomenon cannot be really controlled, most of the time one is led to inject large 
amounts of gas, which requires considerable conipression means. Furthermore, injection 
of large amoimts of gas substantially modifies the gas-oil ratio (GOR), which 
1 0 complicates the phase separation operations at the top of the riser. 

SUMMARY OF THE INVENTION 

The object of the method according to the invention is to exercise, by modelling the 
instability phenomena described above, an eflEicient dynamic control over the pressure 
of the gas to be injected into pipes so as to reduce these phenomena as much as 
15 possible. 

The method according to the invention allows to neutralize, by controlled gas 
injection, the formation of liquid slugs or accumulations at the foot of a pipe portion 
greatly inclined to the horizontal or riser connected to a pipe carrying circulating 
multiphase fluids. This control is essentially exerted by injecting at the base of the riser 
20 a volume of gas substantially proportional to the mass flow rate variation with time of 
the gas phase of the circulating fluids, and preferably substantially equal thereto, when 
this variation is positive. 
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According to another implementation mode, control is exercised by modulating also 
the volume of gas injected by a quantity proportional to the mass flow rate variation of 
the liquid phase of the circulating fluids, also measured with the same time interval. 

Injection is carried out at any time t, from the mass flow rate variation with time of 
5 the gas phase of the circulating fluid, measured at a previous time interval 

The hnplementing device allows to neutralize, by controlled gas injection, the 
formation of liquid slugs or accumulations at the foot of a pipe portion greatly inclined 
to the horizontal or riser connected to a pipe carrying circulating multiphase fluids. It 
comprises gas injection means connected to the base of the riser, means for measuring 
10 the flow rate of the gas phase of the circulating fluids, and a computer designed to 
control the injection, by the injection means, of a volume of gas substantially 
proportional to the flow rate variation with time of the gas phase of the circulating 
fluids, when this variation is positive. 

The con[5)uter is for exan^ple suited to control the injection, by the injection means, 
15 of a volume of gas substantially equal to the flow rate variation with time of the gas 
phase of the circulating fluids. 

According to an embodiment, the device also comprises means for measuring the 
flow rate of the liquid phase circiilating in the pipe, the computer being suited to 
modulate the volume of gas injected by a quantity proportional to the measured flow 
20 rate variation of the liquid phase. 



BRIEF DESCRIPTION OF THE FIGURES 

- Figure 1 diagranmiatically shows a pipeline of low inclination connected to a riser 
greatly inclined to the horizontal, 

- Figure 2 diagramroatically shows a steady flow with continuous penetration of gas 
in a gas-liquid separator at the upper end of a riser, 

- Figure 3 diagrammatically shows the formation of a Hquid slug at the foot of the 
riser, 

- Figure 4 shows the stage when the liquid sli^ reaches the separator at the top of the 
riser, 

- Figure 5 diagrammatically shows the gas pocket getting into the liquid accumulated 
in the riser, 

- Figure 6 shows the stage when the liquid falls back to the base of the riser, 

- Figure 7 shows the pressure curve in the pipeline at the foot of the riser during the 
previous cycle, 

- Figures 8A to 8E diagrammatically and respectively show, in a riser, a flow mode 
referred to as bubble flow (8A), a flow mode referred to as intermittent flow (8B), a 
flow mode referred to as chum flow (8C), a flow mode referred to as annular flow (8D) 
and a flow mode referred to as stranded annular flow (8E), 

- Figures 9A to 9G diagrammatically and respectively show, in a pipeline, a flow 
naode referred to as stratified flow (9 A), a flow mode referred to as wavy stratified flow 
(9B), a flow mode referred to as droplet annular flow {9C), a flow mode referred to as 
dispersed bubble flow (9D), a flow mode referred to as intermittent flow (9E), a flow 
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mode referred to as small-pocket flow (9F), and a flow mode referred to as elongated- 
bubble flow (9G), 

- Figure 10 shows an example of evolution of the pressure at the foot of a 14-m high 
riser, without any control by gas mjection, 

5 - Figures 1 1 A, 1 IB respectively show the flow rate variations of liquid and gas at the 
riser outlet, also without control, 

- Figures 12 and 13A, 13B correspond to Figures 10 and llA, 1 IB for a 250-m high 
riser, 

- Figure 14 shows the evolution of the voltmie fraction of Uquid at the riser outlet, 
10 without control, 

- Figure 15 respectively shows an example of variation of the pressure at the foot of 
the 14-m high riser before gas injection and the stabilization obtained by controlled 
injection according to the method (from 500 s), 

- Figures 16A and 16B respectively show, in the same riser, the curves of the 
15 variation with time of the liquid flow rate and of the gas flow rate, before and after 

control (also from 500 s), 

- Figure 17 shows, in a 250-m high riser, the pressure variation with time at the foot 
of the riser, without control (thick line) and with control (dotted line), 

- Figures 18 A, 18B respectively show, in the same riser, the liquid and gas pressure 
20 variations with time, without control (thick line) and with control (dotted line), 

- Figure 19 shows the compared evolutions of the pressure at the foot of a 14-m riser 
with control of the injection pressure as a function of qo and qL (full line) and only of q© 
(dotted line). 
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- Figures 20A, 20B respectively show, in the same riser, the liquid and gas flow rate 
variations with time, with control of the injection pressure as a fimction of qo and qL 
(dotted line) and only of qo (full line), 

- Figure 21 shows, in a 250-m high riser, the pressure variations with time at the foot 
5 of the riser, with control of the injection pressure as a fimction of q^ and qL (dotted line) 

and only of qo (foil line), 

- Figures 22A, 22B respectively show, in the same riser, the liquid and gas pressure 
variations with time, with control of the injection pressure as a fonction of qo and qL 
(dotted line) and only of qo (foil line), 

10 - Figures 23A, 23B respectively show the respective liquid and gas flow rate 
i|| variations with time at the riser outlet, and 

s| r. 

- Figure 24 diagrammatically shows an embodiment of the device for implementing 

13 

i^j the method, blowing the formation of slugs to be neutralized. 

i3 DETAILED DESCRIPTION 

15 I-l Modelling 

Modelling of the flow phenomena in the system consisting of the pipe and of the 
riser of Figure 1 is obtained by means of the following hypotheses. 

We choose a one-dimensional approximation where we consider the averages of the 
various quantities on a (straight) cross-section of the pipeline or of the riser. Since this 
20 approximation is acceptable only if the radius of curvature of the pipeline is assimied to 
be infinite and its diameter constant, the modelling procedure will concern the parts on 
either side of the connecting elbow. 
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We use a drift flow type modelling procedure with a mass conservation equation 
per phase and wherein the liquid and gas momentum conservation equations are added 
to one another so as to have a single equation referred to as mixture momentum 
conservation equation. To close the system, we choose a relation expressing a friction 
law between the phases. 

We also conventionally assume that : the flow is isothermal, the fluids are 
Newtonian, the gas is a perfect gas, the liquid is incompressible (its density is therefore 
constant) and there is no mass transfer between the two phases. 

We consider that the Mach number of the mixture is small so that the pressure 
waves are propagated at an infmite velocity instead of a velocity close to the sound 
velocity in the mixture. High-frequency phenomena are suppressed but the void fraction 
waves continue to be propagated at a velocity close to the velocity of the mixture. This 
hypothesis is reinforced by the feet that we study the system « in transition » to the state 
of obstruction or disturbed state, i.e. close to the steady state. This hypothesis is 
translated in the model into the absence of inertia terms in the momentum conservation 
equation. 

1-2 Selection of the fundamental quantities 

The following quantities are defined : 

Rg and Rl are the volume fractions of gas and of liquid in the pipes, 
qe and qL are the mass flow rates of gas and of liquid per section unit. 
Vg and Vl are the velocities of the gas and of the liquid. 
pG and pL are the densities of the gas and of the liquid. 
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P is the mean pressure of the mixture. 

All these quantities are fimctions of RxR (space and time), zero on RxR.* and 
continuous on RxR+*. According to the context, it is essential to always precisely say if 
the pipe concerned is the pipeline or the riser. Thus, we work on the pipeline in [0,L]xR 
5 by writing the variables (x,t) and on the riser in [-H,0]3dR by writing the variables (z,t). 

However, if an equation is valid in the pipeline as well as in the riser, the equation 
will be formally written with the variables (x,t) g [-H,L] u R^.. 

1-3 Intrinsic equations 

All these quantities are connected by algebraic and differential relations which do 
10 not depend on the flow considered (these flows are described in the next section), which 
is what we call intrinsic equations. The other types of equation, mainly friction laws, are 
studied in the next section. 

1-3-1 Algebraic equations 

We first express the relations directly obtained from the definition of the quantities. 
15 ^g=/>gVg (I-l) 

(1.3) 

The perfect gas equation allows to establish the following relation between the 
pressure and the density of the gas : 

20 ^^(x,o=^^ a.4) 

a 
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^ or Mg is the molar mass of the gas, T the temperature of the 

mixture and R the perfect gas constant, a corresponds to the soimd velocity in the gas at 
Ibar. 

1-3-2 Differential equations 
5 The naass conservation of each phase imposes that : 

dt dx 
dt dx 

The momentum conservation accordii^ to the previous hypotheses imposes that : 

^ = Sin 0{RoPo + RlPl)- (1-7) 

10 where 9 is the inclination of the pipe, g the gravity constant and Fp the wall friction 
(friction of the stream against the wall). 

A closing equation is added to the aforementioned eqmtions in form of an algebraic 
slippage law as follows : T(P, Rg? Vg, Vl, Pl) = 0. 

1-4 Slippage law 

15 The selected slippage law Y depends on the flow regime. Three flow types can be 

considered : stratified flow, dispersed bubble (or simply bubble) flow and intermittent 
flow. All the flow regimes are illustrated in Figures 8A to 8E and 9A to 9G. For our 
study, we consider the case of an intermittent flow in the riser and of a stratified flow in 
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the pipeline. We shall see hereafter that, considering possible simplifications, no 
slippage law is necessary for the stratified flow regime. 

1-4-1 Intermittent flow 

This flow regime is due to a « superposition » between a bubble flow and a 
5 stratified flow. When the gas flow rate increases, the bubbles clump together and 
coalesce. Large shell-shaped bubbles appear. They are separated by liquid slugs which 
geniKrally contain small gas bubbles. 

Under such conditions, the fidction law for an intermittent flow is expressed as 
follows : 

10 Va -C,{R^V^ ^RjV.yV^ = 0 (LIO) 

hence function ^^nt (Vg, Vl, Rl, V„). 

Besides, Vo^ is experimentally determined and has the following form : 
= (0.35sin 6 + 0.54cos^)VgD 

15 II STUDY OF THE PIPELINE-RISER SYSTEM WITHOUT FRICTION 

In this part, we fix Fp = 0 in Equation 1.7 considering that the riser is vertical and 
that it is assumed that, in the flow disturbances observed, the fi-ictions only have a 
limited influence in relation to the gravity. Solution of the equations is therefore a priori 
simplified. 
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n-l Steady state in the riser 

In the case of a steady state, the various quantities only depend on the space ; they 
are marked with a line to show that they are steady quantities. Furthermore, any 
constant fixnction is asterisked. Thus, for example, under steady conditions, the Uquid 
and gas flow rates and are constant. They are therefore denoted by qL* and qo*- 
Besides, the pressure at the top of the riser is denoted by Pq. 

Intermittent flow 

The closing equation is thus given by 1. 10. Furthermore, as before, in the steady 
state, the gas and liquid flow rates being constant, they are denoted in the same way. 

The implicit formxilation of P under intermittent flow conditions is written as 
follows : 



(11.15) 



with g given by the relation : 



gigc-^plv^+c^qlf 



In 



U'^plK^Co^iy gpl{<ll+PlK^C^qiy 



The pressure can be calculated in implicit form as for a steady flow regime under 
the assumption of an intermittent flow. 



P(z) can thus be assumed to be known at least numerically. It is then posible to 
express all the other quantities as a function of P(z) and z. For the transient states, the 
steady quantities are assumed to be known. 
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n-2 Steady model in the pipeline 

We a priori have the flow regime choice. However, practice shows that it is more 
rational to choose the stratified flow regime, or at a pinch the intermittent flow regime, 
knowing that the second one is less probable. The solution of the intermittent flow 

regime is similar to that of the previous section while not imposing ^ Equation 

1-7. The system can once again be integrated. Assuming that the flow in the pipeline is 
stratified, the closing equation is too oon:5)licated to allow the equations to be solved. 
We therefore involve the pipeline only by the effect it has on the riser, i.e. by measuring 
the riser inlet quantities (flow rates, pressure, surface fi:actions) instead of determining 
them in relation to the calculated pipeline outlet quantities. The calculations are 
therefore replaced by riser inlet measurements. 

II-3 Transient state in the riser in the tangent linear model 

We develop hereafter the calculations about the steady state and linearize the 
intrinsic equations to the first order so as to simply solve the transient system. 

II-3-1 Riemann invariant 

Introduction of this invariant allows to facilitate the pertinent solution of the 



equations in the tangent linear model. It is expressed in the form k = 



which 



has no evident physical significance. We thus have the following relation : 



V(2,0e]-f/,0[x< 



dt 




0 



(n.28) 
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Quantity k moves in the riser, under intermittent flow conditions, at the velocity Vo- 

Thus, the property of k allows us to assert that the propagation of k under 
intermittent conditions is e3q)ressed as follows : 

Kz,t) = k{-H,t-T{z)), (11.30) 

Quantity k at a height z of the riser has the same value as at the bottom of the riser 
T(z), one time unit a earlier. 

With this Riemann invariant, we are going to express P'(-H,t), which is the pressure 
variation at the bottom of the riser where the disturbances start, as a function of the riser 
inflow rates. This can allow us to conceive a control over this pressure. 

II-3-2 Transient infl ow-outflow laws, intermittent flow regime 

We suppose here that the values of the steady quantities are given (these functions 
are known from the steady mode study). By applying the Imearization techniques, we 
can establish the following result : 



P'iz,t) = Fiz) 



f . kC \ 

i^.*UG-H H'l-h (t} + F{-H,t) 

V Pl J 



(n.3i) 



where N, , q'o-H , q'^-H are fimctions of t respectively obtained from N, q'c , q\ wiiere the 
first variable has been set at z, -H and -H respectively. 

In fact, linearization of Equation 1.7 gives : 

= -S\Pg - p1Wg+RgP"o)) (n.32) 
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We are going to replace R'q by its expression as a function of P' and k'. We 
therefore note that k= - By linearizing this relation, we obtain 



^,^^ ^(1 C^Ro) PgRg Thgjjj^ j.gpj^g j^,^ JJ 32 ^ ^ 



obtain the following differential equation : 



ejqjression, we 



The Riemann invariance property of k extends to k' in the following forms 



dk' — dk' — _ 1 

^ + = 0 and k\z,t) = k'(-H,t-T(2)) with T(z) = L^dx 



We can therefore write 



10 Furthermore, we can calculate k' as a fimction of the flow rates, and we obtain ; 



k ^ 



1g 



kC, 



0 



P*L ) 



By injecting this relation into 11.33, this equation becomes 



& a(pc+CoA:) 



SiPG-P*L)(X-CM k 



Pg + 



P*L 



11.34 



«^G 
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We then put : 

a\p^+C^k) Po+Cok go 

This equation is then integrated by means of the constant variation method. We then 
obtain : 



1 



r 



C (x)^q'o i-H,t - Tiz)) - ^g', (-H,t - f{z)) 1 



dx 



J 



Finally, we put 



then 



P\z,t) = F(zi['''M(v) 



g'^(-H,t-v)—^q',i-H,t-v) 



dv + P'i-H,t) 



and if we put N(z,v) = l[(^j.(^)](v)M(v) then : 



10 F{z,t) = Fiz) 

as we have seen above 



* \g'a., '^Z^-" \^ + 



In these relations : 

q'o-H is the mass flow rate variation with time of the gas phase in the circulating 
multiphase fluid at the foot of the riser, i.e. at the height -H in relation to the top of the 
15 riser ; 
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id 



qVn is the mass flow rate variation with time, at the same height, of the liquid phase in 
the multiphase fluid. 

m GAS INJECTION CONTROL 

The action selected to correct the disturbances consists in controlling the pressure at 
5 the bottom of the riser. In fact, if this pressure remains close to a steady value, this 
means that the slugs do not form and that the gas is never really blocked. The action 
will thus concern the pressure P(-H,t) at the bottom of the riser. 

ni-1 Gas lift control 

We assume that the outlet pressure is fixed and therefore that V\ = 0. Relation 
10 111.31 can thus be written as follows : 



No"" 



<1g., r-9i_ 



(0 m-42 



It can immediately be seen that, if one of the two members of the convolution 
product is zero or very small, it is the same for P'(-H,t). Now, No cannot be modified. 



( Ck ^ 
There still is quantity q'G_fj ^q'L^j^ 



\ ~ pL 



(t) , homogeneous with a flow rate, that is 



J 



15 denoted by Q(t). 



The principle of the control mode according to the invention will essentially consist, 
at predetermined action intervals, in measuring at a time ti the above quantity in order to 
obtain a measurement M(ti), then, at the time ti + At, in adding to Q(ti + At) the quantity 
u(ti H- At) = -M(ti), and so on. We thus have Vi > 1 : 
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u{t^ + /A/) = -M{t, + (/ ~ 1) AO 

Q(t^ H- zAO = + m) + + iAt) 

Of course, when At 0, it amounts to equating the second term of the convolution 
product to zero. Since it is physically not conceivable to have a control that sucks the 
gas in a two-phase mixture, we use a control such that : 

u(t^ + iAt) = max(-M(/, + (z - 1)A0, 0) . 
IV VALIDATION OF GAS INJECTION CONTROLS 

To validate this control mode, we use the aforementioned TACITE code which 
simulates flows in the most accurate and realistic way possible. 

We have studied two very different slug formation instances and tested our control 
mode thereon. 

rV-1 Comments on simulations 

The TACITE software uses a finite- volume type method to simulate flows in pipes 
and the pipeline is discretized for example accordmg to the gridding method described 
in patent appUcation FR-EN-00/08,200 filed by the applicant. To simulate gas injection 
at the foot of the riser, we modify the flow between the two grid cells situated just 
before and just after the bend, respectively numbered n - 1 and n. The initial state of our 

C k 

simulations is the steady flow regime. Coefficient is therefore measured once and 

Pl 

for all at the initial time. We identify the steady flow rates with the pipeline inflow 
rates, and our riser inflow rates with the flow rates of grid cell n - 1. To display the 
results, we approximate to the flow of gas to be injected by the flow rate difference 
between grid cells n and n - 1 . 
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IV-2 Simulations of various non-controlled cases 

The simulations carried out by means of the TACITE simulation code are based on 
two geometries where L represents the length of the pipeline, H the height of the riser 
and 0 their common diameter. 

Case 1 :L = 6Om,H=14m,0 = 5cm, 

Case2:L= 1750 m,H = 250m, 0 = 25 cm. 
Case 1 : 14-m riser 

This is the case which we take as the reference : 

Po = 1 bar, Ql = 2.10-2 kg/s and Qq = 2.10-' kg/s. This case has marked oscillations, 
with an oscillation period of the order of one minute. The flow regime corresponds to 
our hypotheses : stratified in the pipeline and intermittent in the riser. 

The oscillations of the system (Figs. 11 A, IIB) are not great enough to completely 
cancel out the flow of gas at the foot of the riser or for the free surfece of the liquid to 
fall below the level of the mouth of the riser. We are therefore not in a situation of 
significant slug formation. This property allows to remam close to our tangent 
linearized hypothesis, while remaining in a flow regime that is not fer from the steady 
state. We therejfore remain, for any t, in the quasi-steady situation of transition towards 
slug formation. 

Two main stages can be distinguished in the cycle : 

1 . A stage of accumulation of the Uquid in the riser, the liquid flow rate is zero at the 
outlet and the pressure increases ; 
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2. A stage of liquid slug production where the pressure decreases. 
Case 2 : 250-m riser 

The main characteristics of this case are as follows : 
Po = 1 0 bar, Ql 4 kg/s and Qg = 0.5 kg/s. 

These conditions are closer to the real operating conditions of a pipeline. In the 
absence of any control, the system enters into a liquid slug formation stage. It can be 
observed that the liquid flow rate at the outlet goes through already violent expulsion 
stages, even before the slug expulsion and pressure fell stage. At the end of the 
simulation, the pressure has reached a maximum value, and the system is about to enter 
into the liquid slug expulsion stage. The system is therefore in the phase of transition to 
slug formation only during the first moments, because the starting point of the 
simulation is the steady state. 

Figure 14 shows very fast cycles in the evolution of the liquid fraction. This is a 
sign of high instability in this case. 

IV-3 Control of both cases 

The first control that we are going to study is the theoretical control found in the 
previous chapter. Following our observations, we shall see that it is possible to conceive 
another control, that we are also going to test. 

IV-3-1 Theoretical control test 

We introduce here a control by gas lift by injecting, at the foot of the riser, a flow of 
gas of the following form : 
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«(0 = max(a(^'^^-^^* )(t-At)) 

P L 

Case 1 : 14-in riser 

We start control just before t - 500 s. The graphs of Figs.IV-16, IV-17 show in 
parallel the evolution of the free system (thick hne) and that of the controlled system 
5 (fine dotted line). This control allows to maintain the pressure and the outlet flow rates 
close to their steady value. The mass of gas injected is small in relation to what comes 
from the pipeline ; the outlet flow rate increase due to the injection of gas is less than 
5%. 

Case 2 : 250-m riser 

10 For this simulation, we start control from t = 0 s. The graphs of Figs.15 and 16A, 

16B show in parallel the evolution of the free system (thick line) and that of the 
controlled system (fine dotted line). The gas flow rates in the presence of control are 
much more regular than in the previous case. 

Control by gas injection using the theoretical formula thus functions in both cases, 
15 although they are very different. 

IV-3'2 Simohfied control 

Since a great correlation is observed between the evolutions of the gas and liquid 
flow rates, and since we almost always have j^*^^ ^ | > 
control the system with the simplified control : 
20 u{t) = max(0, {q'^ ^ )(t ^ At)) 



, it is also possible to 
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Case 1 : 14-m riser 

In the graphs of Figs. 17, 18 A, 18B, control with qo and is shown with a thick 
line, whereas control with qo only is shown in fine dotted line. Control starts from t = 
Os. 

5 The results of the control concerning only qo are practically identical to those 

obtained with the control already tested, or even slightly better as regards the pressure 
oscillations and the control speed. 

Case 2 : 250-m riser 

In this distinctly more imstable case, we observe (Figs.22A, 22B) no fest 
10 oscillations of the gas flow rate with the simplified control The system controlled with 
qo evolves with a slight lead in relation to the same case controlled with the conqjlete 
e?q)ression. 

We can therefore see that the simplified control without coefficient qL also allows to 
control the system, but the theoretical control is better because, unlike the simplified 
15 control, it does not lead to problems at the riser outlet, it is more economical as regards 
injection gas and the system is controlled just as well. 

IV-3-3 Control robustness 

We test the robustness of our theoretical control (qo and qO in relation to the 
reaction time of the sensors and of the actuators. We therefore compare the control 
20 obtained above by adjusting the gas flow rate at each time interval (thick dotted line) 
with the control obtained by adjusting these parameters with a lower frequency (fine 
dotted line). 
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Case 1 : 14-m riser 

Control is readjusted every 3 seconds here, which corresponds to 20 time intervals 
for TACITE. Control starts from t = 0 s. 

It can be seen in Figs.23A, 23B that the flow rates are a little more irregular, but the 
5 system is still controlled. 

Case 2 : 250-m riser 

The system is not stabilized with a 2-second time constant. With a 1 -second tnne 
constant, the system is stabilized but we remain close to the order of magnitude of the 
calculating interval (about 0.3 s in this case). Furthermore, there is practically no 
10 difference with the case where control is adjusted at each time interval 

In any case, we have a good control even though the flow rates oscillate very 
slightly around a mean value. 

The device for implementing the method con^rises (Fig.24) gas injection means 1 
connected to the base of the riser, means 2 for measuring the flow rate of the gas phase 
15 of the circulating fluids, and a computer 3 intended to control injection, through 
injection means 1, of a volume of gas proportional to and preferably substantially close 
to the flow rate variation with time of the gas phase of the circulating fluids, when this 
variation is positive. 



